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Abstract 

We introduce a novel multiobjective optimization algorithm based on the 
conformational space annealing (CSA) algorithm, MOCSA. It has three 
Q-i characteristic features: (a) Dominance relationship and distance between 

solutions in the objective space are used as the fitness measure, (b) update 
rules are based on the fitness as well as the distance between solutions in the 
decision space and (c) it uses a constrained local minimizer. We have tested 
MOCSA on 12 test problems, consisting of ZDT and DTLZ test suites. 
Benchmark results show that solutions obtained by MOCSA are closer to 
the Pareto front and covers a wider range of the objective space than those 
r-i ■ by the elitist non-dominated sorting genetic system (NSGA2). 

P-i 

Keywords: conformational space annealing, multiobjective optimization, 
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1. Introduction 

The multiobjective optimization problem (MOOP) is to optimize two or 
more objective functions simultaneously, subject to given constraints. The 
multiobjective optimization can be applied to problems where the final deci- 
sion should be made considering two or more conflicting objectives. MOOP 
occurs in various fields such as industrial design, finance, management and 
many engineering areas. Practical goals in these fields can be generalized in 
such a way that the cost of a process is minimized while the quality of its 
product is maximized. The primary goal is to find a set of solutions that 
any individual objective function cannot be improved without deteriorating 
the other objective functions, and such a set is called a Pareto set. For 
efficient decision making, a set of generated solutions (QS) should meet two 
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conditions: It should be as close to the Pareto front as possible and the 
solutions should be distributed as widely as possible. 

Evolutionary algorithm ( EA) is one of the most popular and successful 
approaches to solve MOOPs dDebL l200ll : fCoello Coello et all 120071 ) . A num- 
ber of EA-based algorithms have been suggested including the vector evalu- 



ated genetic algorith m (VEGA) f Schafferl . Il985l ). the niched Pareto genetic 



algorithm (NPGA) flHorn et al. 
algorithm II fNSGA2) flDeb et al 



1994!) the nondominated sorting genetic 
' l2000h. the strength Pareto evolutionary 



algorithm II (SPEA2) (IZitzler et all 120011). the mimet ic Pareto archived 



evolution strategy fM-P AES) flKnowles and Cornel. l2000l) and mic ro genetic 
algorithm (micro-GA) (ICoello Coello and Toscano Pulidol . 120011 ) . Among 
them, NSGA2 and SPEA2 are arguably the most w idely used methods - 



Other approaches in c lude simulated ann ealing (SA) (ISuman and Kumar 



2005 



Nam and Park . 2000 ). tabu search ( Hansen . 1997:~lGandibleux et al 



1997 D, particle swarm optimizatio n (PSO) flParsopoulos and V rahatis. 



2002 



Coello Coello and Lech uga. 20 Q2|), immune alg orithm (IA) (ICoello Coello and Cortes 



2002 



2004 



20031 ) 



Gao and W ang. 



2011 



Luh et all 120031 ). ant syst em (jDoerner et al. 



Baran and Schaererl . 120031 ) and cultural algorithm (ICoello Coello and Becerral . 



Conformational Space Annealing (CSA) is a highly efficient single-objective 
global optimization algorithm which incorporates advantages of genetic al- 
gorithm and SA. It has been successfully applied to diverse single-objective 



rotein structure 



(Lee et al.. 


1999: 


Pillardv et al.. 


2001; 


Liwo et all 


1999; 


Joo et al. 



2009), finding the minimum energy so lution of a Lenard - Jones cluster fjLee et al. 



20031 ). multiple seq uence alignment f lJoo et all 120081 ) and the community 



detection problem (ILee et all 120121 ) on networks. In these studies, CSA is 
shown to perform more efficient sampling using less computational resources 
than the conventional Monte-Carlo (MC) and SA methods. 

Here, we introduce a new multiobjective optimization algorithm by using 
CSA, MOCSA. Compared to existing EAs, MOCSA has the following dis- 
tinct features: (a) The ranking system considers the dominance relationship 
and the distance between solutions in the objective space, (b) solutions are 
updated by using a dynamically varying distance cutoff measure to control 
the diversity of the sampling in the decision space, and (c) a gradient-based 
constrained minimizer is utilized for local search. 

The remainder of this paper is organized as follows. In section 2, the 
definition of MOOP and related terms are described. In section 3, details 
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of MOCSA is presented. Numerical results and the comparison between 
MOCSA and NSGA2 on various test problems are presented in section 4. 
The final section contains the conclusion. 

2. Problem statement 

The mathematical definition of a MOOP can be defined as follows, 
min v = min f (u) = min (fi (u) , f 2 (u) , 

u u u 

s.t. u = (m, U2, ■ ■ ■ , u n ) G U C R n 

v = (ui,u 2 ,-..,um) eV cR m 

where u is the decision vector, U the decision space, v the objective vector 
and V the objective space. Due to the presence of multiple objective func- 
tions, a final solution of MOOP consists of a set of non-dominated solutions 
instead of a single point. The notion of dominance and related terms are 
defined below. 

Definition 1. A decision vector Ui is said to dominate another solution u 2 
(denoted by ui -< u 2 ), if and only if 

Vi G {1, . . . , m} : /,( Ul ) < /,(u 2 ) A3k G {1, . . . , m} : / fc ( Ul ) < / fc (u 2 ). (1) 

Definition 2. A solution u is said to be non-dominated by any other so- 
lutions (a Pareto optimal solution) if and only if 

^3u* G U : u* ■< u. (2) 

Definition 3. For a given MOOP, a Pareto optimal set in the decision 
space, VS, is defined as 

VS := {u|-3u* Gf/:u^u}. (3) 

Definition 4. For a given MOOP, a Pareto optimal set in the objective 
space, VJ 7 , is defined as 

VT:= {f(u)\ueVS}. (4) 

Since the size of Pareto optimal front, VT is infinite in general, which 
is impossible to obtain in practice, practical algorithms for MOOP yield a 
set of non-dominated solutions of a finite size. It should be noted that VJ 7 
is always a non-dominated set by definition while a non-dominated set of 
solutions generated by an algorithm, which is denoted as a QS, may not be 
a subset of VJ '. 
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...,/ m (u)), 



3. Description of conformational space annealing 

Here, a new multiobjective optimization algorithm based on CSA is de- 
scribed. The CSA was initially developed to obtain the protein structure 
with the minimum potential energy, i.e., to solve a single objective opti- 
mization problem. CSA has been successfully applied to various kinds of 
optimization problems with modification. The general framework of CSA is 
shown in Figured! and the description of MOCSA is given in Algorithm [TJ 



Initial N random 
solutions 



Constrained 
minimization 



Add to First Bank 



Dcut = Dave/2 




Figure 1: Flow chart of CSA is shown. 
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Algorithm 1 Multiobjective CSA algorithm 

1: procedure MOCSA(iV, JV fl , JV C , N m , G = 300, a) 

2: Initialize the bank, P, with N random individuals 

3: Minimize(P) using a constrained local minimizer 

4: Initialize seed flags of all individuals to zeros : s(i) 0, % — 1, . . . , N 

5: Get average distance, D ave , between all pairs of individuals and set 
D cut as D ave /2: D cut <— D ave /2 

6: Initialize generation counter to zero : g 

7: Initialize the reserve bank, R, to an empty set 

8: while g < G do 

9: if s(i) == 1 Mi then 

10: Generate random individuals, P' 

11: Minimize (P') 

12: P^PHP' > Expand search space 

13: end if 

14: Evaluate Fitness of P 

15: Select N s seeds among individuals with s(i) = and set s(i) to 

1 

16: T c ^— Generate iV s A^ c trial solutions by crossover 

17: T rn ^— Generate N s N m trial solutions by mutation 

18: T^T c UT m > Trial solutions 

19: if G%5 == then 

20: Minimize(T) 

21: end if 

22: Update(P, T, M) 

23: D cut i- max(aD cut ,D ave /5) > Reduce D cut 

24: g <- g + 1 

25: end while 
26: end procedure 
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3.1. Conformational space annealing 

CSA is a global optimization method which combines esse ntial ingredi- 
ents o f three methods: Monte Carlo with minim ization f MCM) flLi and Scheraeal . 



1987h . genetic algorithm (GA) flGoldbereJ . [l989h . and SA flKirkpatrick et al. 



19831 ). As in MCM, we consider only the solution/conformational space of 
local minima; in general, all solutions are minimized by a local minimizer. 
As in GA, we use a set of N solutions (called bank in CSA, denoted as P) 
collectively, and we generate offsprings from the bank solutions by cross- 
over and mutation. Finally, as in SA, we introduce a distance parameter 
D cut , which plays the role of the temperature in SA. In CSA, each solution is 
assumed to represent a hyper-sphere of radius D in the decision space. Di- 
versity of sampling is directly controlled by introducing a distance measure 
between two solutions and comparing it with D cut , to prevent two solutions 
from approaching too close to each other in the decision space. Similar to 
the manipulation of temperature in SA, the value of D cut is initially set to 
a large value and is slowly reduced to a smaller value in CSA; hence the 
name conformational space annealing. 

Compared to the conventional EA for multiobjective problems, MOCSA 
has three distinct features; (a) a ranking algorithm which considers the dom- 
inance relationship as well as the distance between solutions in the objective 
space, (b) an update rule with a dynamically varying distance cutoff mea- 
sure to control the size of search space and to keep the diversity of sampling 
in the decision space and (c) the usage of a gradient-based constrained min- 
imizer, feasible sequential quadratic programming (FSQP), for local search. 

In CSA, we first initialize the bank, P, with N = 50 random solutions 
which are subsequently minimized by FSQP constrained minimizer. The 
solutions in the bank are updated using subsequent solutions found during 
the course of optimization. The initial value of D cut is set as D avg /2, where 
D avg is the average distance in the decision space between two solutions at 
the initial stage. A number of solutions (20 in this study) in the bank are se- 
lected as seeds. For each seed, 30 trial solutions are generated by cross-over 
between the seed and randomly chosen solutions from the bank. Additional 
5 are generated by mutation of the seed. It should be noted that if a solu- 
tion is used as a seed and not replaced by a offspring, it is excluded from the 
subsequent seed selection. The generated offsprings are locally minimized 
by FSQP which guarantees to improve a subset of objective functions with- 
out deteriorating the others and without violating given constraints. To 
limit the computational usage, the minimization is performed only once per 
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every five generation steps (see Algorithm [T]) . 

Offsprings are used to update the bank, and detailed description on the 
updating rule is provided in Section 13.31 Once all solutions in the bank 
are used as seeds without generating better solutions, implying that the 
procedure might have reached a deadlock, we reset all bank solutions to be 
eligible for seeds again and repeat another round of the search procedure. 
After this additional search reaches a deadlock again, we expand our search 
space by adding additional 50 randomly generated and minimized solutions 
to the bank (N = N + 50), and repeat the whole procedure until a termi- 
nation criterion is satisfied. The maximum number of generation is set to 
G. Typically, with G = 300 in this study, MOCSA is terminated before a 
deadlock occurs with the final bank size of N = 50. 

3.2. Fitness function 

For a given set of generated solutions, QS, the fitness of solution i is 
evaluated in terms of rij, and d\ 2 . rij is the number of solutions in QS 
which dominate i. reij is the number of solutions in QS dominated by i. d} 2 
is the sum of distances from i to its nearest and second nearest neighbors in 
QS in the objective space. The relative fitness between two solutions, i and 
j, is determined by the comparing function shown in Algorithm |2j With 
a set of non-dominated solutions all values of rij and m; become zeros and 
the solution with the least value of d 12 is considered as the worst. 

3.3. Update rule 

The solutions generated by crossover and mutation are locally minimized 
by FSQP constrained minimizer and we call them trial solutions. Each trial 
solution t, is compared with the bank P for update procedure as shown 
in Algorithm [3j First, w, the closest solution in P from t in the decision 
space is identified. If there exist dominated solutions in P, the closest con- 
formation search is performed only among them. Otherwise, P is a set of 
non-dominated solutions, and all in P are considered. Once w is found, 
the distance D in the decision space between t and w is calculated. If 
D > D cut , the current cutoff distance, which indicates that t lies in a newly 
sampled region in the decision space, remote from the existing solutions in 
P, the dominance relationship between t and the worst solution in P, u, is 
compared. If D < D cut , t is compared with w. The selection procedure, 
described in Section I3.4[ is performed to determine which solution should 
be kept in P. At each iteration step, D cut is reduced with a pre-determined 
ratio, a. After D cut reaches to its final value, D avg /5, it is kept constant. 
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3.4- Selection procedure 

In Algorithm El for a given trial solution t and a solution x in P, P is 
updated as follows. If t dominates x, t replaces x. If t is dominated by x, t 
is discarded and x stays in P. If there is no dominance relationship between 
t and x and if x is better than t by Algorithm [2j t is discarded and x 
stays in P. Finally, when t is better than x without dominance relationship 
between them, Algorithm @] is used. 

For the selection procedure, we introduce an additional set of non- 
dominated solutions, the reserve bank, R. Due to the limited size of the 
bank, we may encounter a situation where a solution exists in P which is 
not dominated by the current bank, but dominated by a solution elimi- 
nated from the bank in an earlier generation. To solve this problem, non- 
dominated solutions eliminated from P are store d up to 500 in R, which is 
conceptually simila r to an archive in other EAs (IKnowles and Cornel . I200CH ; 
Zitzler et all 120011 ) . The difference is that K is used only when more than 
half of the solutions in P are non-dominated solutions because CSA focuses 
more on diverse sampling rather than optimization at the early stage of the 
optimization. Note that P keeps only non-dominated solutions. 



4. Test suites 

For the benchmark test of MOCSA, we have selected 12 widely used 
test pr oblems in the fiel d. They consist of ZDT ( iZitzler et all 120001 ) and 
DTLZ (IDeb et all 120021 ) . Each test suite contains several functional forms 
and can feature various aspects of optimization algorithms. Comprehensive 
analysis on the characteristics of the two test suites are well documented by 
Huband et al. fjHuband et al.l . 120061 ) . In both suites, the input vector, x, is 
divided into two sets y and z to construct test problems as follows, 

Given x = {xi, . . . , x n } 

let y = {yi,...,y,} = {x l ,...,x j } 
z = j>i, ...,z k } = {x j+1 , . . . ,x n } 

, where n and j are the dimensions of decision and objective spaces respec- 
tively and k = n — j. 



4.1. ZDT 

The ZDT problem suite consists of six test problems and is probably the 
most popular test suite to access multiobjective optimization algorithms. 



Algorithm 2 Comparing function 
1: function COMPARE^, j) 

2: 7ij(j) Number of solutions in P which dominate i(j) 

3: m i(j) Number of solutions in P which are dominated by i(j) 

4: d\?--s The sum of distances from i to the nearest and second nearest 

neighbors in P in the objective space 
5: if rii < rij then 
6: return i is better 

7: else if rii > rij then 
8: return j is better 

9: else > rii == n j 

10: if rrii > rrij then 

11: return i is better 

12: else if m, < rrij then 

13: return j is better 

14: else > mi == rrij 

15: if df > df then 

16: return i is better 

17: else 

18: return j is better 

19: end if 

20: end if 

21: end if 

22: end function 
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Algorithm 3 Update procedure 



1: procedure Update(P, T, R) 

2: for t in T do 

3: if all solutions in P are non-dominated then 

4: w ^— Nearest solution to t in P in the decision space 

5: else 

6: w ^— Nearest solution to t among dominated solutions in P 

in the decision space 

7: end if 

8: if distance(t, w) > D cut then > distance in the decision space 

9: u Worst solution in P 

10: x <- u 

11: else 

12: X <r~ W 

13: end if 

14: if t -< x then > t dominates x 

15: t replaces x 

16: else if t y x then > x dominates t 

17: x stays in P and t is discarded 

18: else if x is better than t by Algorithm [2] then 

19: x stays in P and t is discarded 

20: else 

21: SELECT(t,x) 

22: end if 

23: end for 

24: end procedure 
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Algorithm 4 Select procedure 



1: procedure SELECT~(t,x) > t G T, x G P, t is better than x 

2: if \J\fV(F)\ < |P|/2 or x is dominated by x' G P then 
3: t replaces x in P 

4: else > x is non-dominated in P and R is used 

5: if t is not dominated by R then 

6: t replaces x in P 

7: else > t is dominated by R 

8: u ^— Nearest dominating solution to t in R in the objective 

space 

9: Move u from R to P 

10: end if 

11: Move x from P to R 

12: R <- W£>(R) 

13: end if 
14: end procedure 



The explicit functional forms of five ZDT problems are presented in Ta- 
ble [TJ The ZDT test suite has two main advantages: (a) The Pareto fronts 
of the problems are known in exact forms and (b) benchmark results of 
many existing studies are available. However, there are shortcomings: (a) 
The problems have only two objectives, (b) none of the problems contain 
flat re gions and (c) none o f the problems have degenerate Pareto optimal 



front flHuband et all 120061 ) 



4.2. DTLZ 

The original DTLZ suite consists of nine test problems which are scalable 
to any number of objectives. This scalability is important since it makes the 
test suite suitable for testing algorithms for many objective problems. The 
explicit functional forms of the first seven DTLZ problems are presented in 
Table [2J DTLZ8 and DTLZ9 are omitted in many benchmark studies due 
to their additional constraints and they are also omitted in this study. 



5. Performance measures 

To evaluate the performance of an algorithm, two aspects are considered: 
Solutions should be (a) as close to the Pareto front as possible and (b) as 
diversely distributed among them. Here, we have used four measures defined 
below. 
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Table 1: Five real- valued ZDT problems are described. The first objective depends only 
on the first decision variable as /i(yi) and the second objective is given as f '2(2/1, z) = 
#( z )M/i(2/i),3( z )), where y = {y u ...,y j } = {x!,...,Xj} and z = {zi,...,z k } = 
{xj + \, . . . ,x n }, where n and j are the dimensions of the decision space and the ob- 
jective space. Unless the functional forms of f±,g and h are separately given, they are 
identical to those of ZDTl. 



Name 


Problem 


domains 


ZDTl 


h = Vi 

■ k 

< 'i=\ 

Kfi,g) = 1 - Vfi/g 


[0,1] 


ZDT2 


h = l- {h/gf 


[0,1] 


ZDT3 


h = l- y/fjg - (f l /g)stn(10nf 1 ) 


[0,1] 


ZDT4 


g = l + 10k + J2ti( z i - 10cos(4irzi)) 


yi e [0, 1] 

Zi, . . . , z k G [—5, 5] 




/x = 1 - exp(-4i/i)sm 6 (67ri/i) 




ZDT6 


k 

g = l + 9(Y,^/k)) - 25 
h = l- {h/gf 


[0,1] 
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Table 2: Seven real-valued DTLZ test problems are described. The objective space is 
j-dimcnsional and j is set to 3 for benchmarking in this study. The input vectors are 
y = {yi,...,yj} = {xi,...,xj} and z = {zi,...,z k } = {x j+ i, . . . ,x n }, where n is the 
dimension of decision space and k — n — j. 



Name 


Problem 


domains 


DTLZ1 


1 ' _1 

i=i 

1 /j-m \ 
fm=2:j-l = -(l + g) \ T[yij (1- Vj-m+l) 

f j = l(l + g)(l- yi ) 

k 

g = 100[k + ((*• ~ 0-5) 2 - cos(20n(zi - 0.5)))] 
i=i 


[0,1] 


DTLZ2 


1—1 

A = (1 + c) TT COs(tT%/2) 

1=1 

/m=2:j-i = (1 + fiO I JJ cos(7ryi/2) J sm(7ry j _ m+ i/2) 
/j = (! + g)sin('Ky 1 /2) 

k 

g = - 0.5) 2 
i=i 


[0,1] 


DTLZ3 


Same as DTLZ2, except (7 is replaced by the one from DTLZ1 


[0,1] 


DTLZ4 


Same as DTLZ2, except yi are replaced by y° -5 


[0,1] 


DTLZ5 


Same as DTLZ2, except 2/2,. ..j-i are replaced by 2(1+3) ' 


[0,1] 


DTLZ6 


Same as DTLZ5, except (7 is replaced by (7 = X^Li ^i' 1 - 


[0,1] 


DTLZ7 


fm=l:j—l 2/m 

£ = (! + </) " E[j^(l + «n(37r/i))l j 
i=i 


[0,1] 



5.1. Spacing (S) and average distance between solutions ((d)) 

The spacing, S, measures how uniformly generated solutions QS are 
distributed in the objective space and it is defined as 



S 



where di 



i 



N 



/ _ fk 



i=l 
m 



min > I f. 

N 



(d) = j2 d i/ N 



i=l 



, where N is the number of solutions in QS, di is the distance between the 
solution % and its nearest neighbor in the object ive space in terms of the 
Hamming distance and (d) is the average of di (jScottl . Il995l ). Ideally, we 
want a larger value of (d) with a smaller value of S. 

5.2. Generational distance (GD) 

The generational distance measures the average distance between the 



gener ated solutions {QS) and the Pareto front (IVan Veldhuizen and Lamont 



2000|). GD is defined as 



GD 




(5) 



, where N is the number of solutions in QS and U is the Euclidean distance 
between the solution i £ QS and its nearest in VT '. If all solutions in QS 
lie on VJ 7 , GD is 0. Therefore, a lower GD value is preferred. 

5. 3. Error ratio ( ER ) 

The error ratio measures the fraction of solutions which are not on VT 
and it is defined as 



N 



ER = - (6) 

, where iV is the number of generated solutio ns and Cj = if the solution i 
belongs to VT within 0.01, = 1 otherwise (IVan Veldhuizen and Lamont 
l2000h . 
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Table 3: Dimensions of decision and objective space of benchmark problems are shown. 



Problem 


dimension of decision space 


dimension of objective space 


ZDT1-3 


30 


2 


ZDT4,6 


10 


2 


DTLZ1-7 


10 


3 



6. Result 

The numerical results on 12 test problems obtained by MOCSA and 
NSGA2 are shown in Table H] and Figures [2] & [31 For benchmarking test, 
the number of population and trial solutions were set to 50 and 700 for both 
methods. The NSGA2 was iterated for 250,000 generations and probabil- 
ities of crossover and mutation were fixed at 0.9 and 1/24. For MOCSA, 
the number of seeds used to generate trial solutions was set to 20. The di- 
mensions of decision and objective space of benchmark problems are listed 
in Table [3J 

For ZDT problems, MOCSA outperforms NSGA2 in terms of both con- 
vergence and diversity of the solutions, except ZDT6. For ZDT6, MOCSA 
is better than NSGA2 by (d) but worse by S. For all ZDT problems, entire 
solutions obtained by MOCSA lies within 0.01 on the Pareto front, shown 
in red circles in Figure [21 while NSGA2 solutions are slightly off as shown 
in Table [U Regarding the spread of the solutions, MOCSA results are more 
uniformly distributed on the Pareto front, signified by lower S values, than 
NSGA2 for the first four ZDT problems. 

For DTLZ problems, MOCSA outperforms NSGA2 in almost all aspects. 
For DTLZ7, 8% of solutions by MOCSA is not on the Pareto front, while 
22% is the case for NSGA2. In terms of S and (d), MOCSA provides more 
evenly distributed solutions, covering a wider range of the objective space. 
It should be noted that (d) by MOCSA is larger than that by NSGA2 for for 
all 12 problems tested (by the factor of 1.59 on average), which is signified 
by the fact that red circles by MOCSA cover Pareto fronts evenly while blue 
crosses by NSGA2 appear to be locally and unevenly clustered in Figures [2] 

7. Conclusion 

In this paper, we have introduced a novel multiobjective optimization 
algorithm by using the conformational space annealing (CSA) algorithm, 
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Table 4: Four performance measures are shown for MOCSA and NSGA2 



Problem 


MOCSA 


NSGA2 


(d) 


S 


GD 


ER 


(d) 


S 


GD 


ER 


ZDT1 


0.0404 


0.0055 


0.0000 





0.0270 


0.0156 


0.0011 


0.04 


ZDT2 


0.0404 


0.0082 


0.0000 





0.0292 


0.0146 


0.0212 


0.02 


ZDT3 


0.0438 


0.0148 


0.0001 





0.0329 


0.0201 


0.0020 


0.02 


ZDT4 


0.0404 


0.0097 


0.0000 





0.0328 


0.0159 


0.0006 


0.02 


ZDT6 


0.0327 


0.0150 


0.0000 





0.0216 


0.0119 


0.0000 





DTLZ1 


0.1114 


0.0068 


0.0000 





0.0615 


0.0319 


0.0000 





DTLZ2 


0.2319 


0.0646 


0.0021 


0.02 


0.1361 


0.0683 


0.0020 


0.04 


DTLZ3 


0.2770 


0.0225 


0.0000 





0.1139 


0.0739 


0.0000 





DTLZ4 


0.2478 


0.0424 


0.0009 





0.1630 


0.0898 


0.0019 


0.02 


DTLZ5 


0.0487 


0.0059 


0.0000 





0.0309 


0.0176 


0.0610 


0.06 


DTLZ6 


0.0484 


0.0156 


0.0000 





0.0306 


0.0135 


0.0000 





DTLZ7 


0.2897 


0.0510 


0.0011 


0.04 


0.1880 


0.1322 


0.0071 


0.22 



MOCSA. Benchmark results on 12 test problems show that MOCSA finds 
better solutions than NSGA2, in terms of four criteria tested. Solutions by 
MOCSA are closer to the Pareto front, a higher fraction of them are on the 
Pareto front, they cover a wider objective space, and they are more evenly 
distributed on average. We note that the efficiency of MOCSA arises from 
the fact that it controls the diversity of solutions in the decision space as 
well as in the objective space. 
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